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The glassy behavior observed in the pyrochlore magnet Y2M02O7, where the magnetic Mo 4+ 
ions interact predominantly via isotropic nearest neighbor antiferromagnetic exchange, possibly 
with additional weak disorder, is a distinct class of spin glass systems where frustration is mostly 
geometrical. A model proposed to describe such a spin glass behavior is the Heisenberg model on 
a pyrochlore lattice with random but strictly antiferromagnetic exchange disorder. In this paper, 
we provide compelling numerical evidence from extensive Monte Carlo simulations which show that 
the model exhibits a finite temperature spin glass transition and thus is a realization of a spin glass 
induced by random weak disorder from spin liquid. From our results, we are led to suggest that the 
spin glass state of Y2M02O7 is driven by effective strong disorder. 

PACS numbers: 75.40.Cx,75.40.Mg,75.50.Lk 



Most magnetic materials develop long range magnetic 
order when the temperature is sufficiently low. There 
are, however, two prominent exceptions: spin liquid (SL) 
and spin glass (SG) systems. Both of them are commonly 
found among frustrated systems where the ordering ten- 
dency is reduced. The SL usually occurs in geometrically 
frustrated systems, whereas the SG arises in random frus- 
trated systems due to competing random antiferromag- 
netic (AFM) and ferromagnetic (FM) couplings^. 

The Y2M02O7 pyrochlore Heisenberg antifcrromag- 
net2r^, possibly with some form of weak disorder, whose 
origin remains both mysterious and controversial^—, does 
not fall in the category of conventional SG materials^. 
The Heisenberg AFM model on the three-dimensional 
pyrochlore lattice of corner sharing tctrahedra is a clas- 
sical spin liquid (CSL) with macroscopically degenerate 
ground states which satisfy the zero net magnetic mo- 
ment (<J>t = 0) constraint on each tetrahedron^. This con- 
straint leads to a gauge theory description of the CSL and 
the prediction of a power-law decaying spin-spin correla- 
tion function of "dipolar" form£. Thus, a clean Heisen- 
berg AFM model on the pyrochlore lattice does not dis- 
play a SG state. However, the macroscopic degeneracy in 
this model opens up the interesting possibility that weak 
random disorder in the spin-spin exchange interactions, 
so small that no competing AFM-FM coupling is present 
in the bare Hamiltonian, may be sufficient to induce a SG 
phase at nonzero temperature. The main question that 
we address in this paper is whether dense random weak 
disorder in the AFM exchange can induce a thermody- 
namic SG phase from a CSL as a case study of what may 
be occurring in Y2Mo20?2~— . 

We study the Hamiltonian H defined on a pyrochlore 
lattice first proposed by Bellier-Castella et aL— : 



rl — Ho + Tidi, 



(1) 



with Ho = Jo E(i,i> Si ■ Sj and H dis = J% s S t ■ 

Sj, where {Si} are three component unit vectors and 



summations are over nearest neighbors, rlo is the usual 
pyrochlore Heisenberg AFM model which displays on its 
own a CSL£ while Hdis describes the random disorder 
which mimics the situation in Y2Mo20^. Weak disorder 
here means | J^ ls | <C Jo- We set the Boltzmann constant 
fcs = 1, and also set Jo = 1 which serves as overall 
energy scale. The bond disorder is uniformly distributed 
in the range jff = [-W, W] with W = J o /10 used in 
the calculations. We refer to H in Eq.(l) as the BGHM 
models. 

The original work of Bellier-Castella et al. suggested, 
on the basis of measurements of the SG overlap parame- 
ter—, that glassy behavior in TL sets in at a temperature 
roughly the same as W— . In more recent studies, Saun- 
ders and Chalker— and Andreanov et aiw^ computed the 
SG correlation function and SG susceptibility. Based on 
numerical data and analytic arguments, the authors of 
Ref. [l2| suggested that there exists a thermodynamic SG 
transition at a nonzero temperature for arbitrary small 
but nonzero W. 

SG simulations techniques have significantly im- 
proved over the past ten years or soi^r— . Recent 
extensive Monte Carlo simulations employing these 
improvements^— have led to the revision of the old be- 
lief that the lower critical dimension, di , of the Edwards- 
Anderson (EA) isotropic Heisenberg SG model is above 
threei. There is no rigorous analytic approach to deter- 
mine the di of the BGHM modeller— and the previous 
MC simulations^— do not come close to the computa- 
tional standard of recent studies of the EA Heisenberg 
SG^£~— . Thus, the analytic arguments and numerical 
data at hand can hardly provide convincing evidence for a 
thermodynamic SG phase in T~L. The BGHM model, with 
its underlying CSL state in the disorder-free regime, as 
well as its broad relevance to the SG behavior observed 
in numerous geometrically frustrated magnetic materi- 
als, make it a model of fundamental significance in the 
field of frustrated magnetism. It is therefore important 
to carefully assess whether rl sustains a thermodynamic 
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SG phase at nonzero temperature, and to reach such a 
conclusion on the basis of numerics that approach the 
standard of SG simulations of EA models^—. 

We first summarize the details of our MC simula- 
tions. A 16 site cubic unit cell for the pyrochlore lat- 
tice is used for generating cubic simulation cells with 
N = 16L 3 spins with L = 4, 6, 8. One Metropolis 
sweep, 2L over-relaxation sweeps and one parallel tem- 
pering swapi^ is defined as one elementary MC step 
(MCS). The temperatures explored for each simulation 
are T^ n > = T min a n where T min is the lowest temperature 
considered and n £ [Q,Nt — 1], where Nt is the num- 
ber of thermal replicas. Thus the highest temperature is 
T max = Tmino^- 1 and a = \/T max /T min . The er- 
ror bars are sample-to-sample fluctuations calculated via 
the jackknife method. Table fl] lists the parameters used 
in our MC simulations. 



L 




Tmax 


N T 


srcquil. 

JV MCS 


a j measurement 
JV MCS 


^ samples 


4 


0.012 


0.028 


24 


2 x 10 4 


2 x 10 4 


1750 


6 


0.012 


0.028 


30 


4 x 10 4 


2 x 10 4 


1763 


8 


0.012 


0.028 


46 


8 x 10 4 


2 x 10 4 


1613 



TABLE I: Parameters of the Monte Carlo simulations. 



To characterize a putative SG phase, we use a param- 
eter defined as the overlap between two thermal replicas 
with the same realization of random couplings, {Jfj S }, 

1 N 

^(k) = S S S S ex P(* k • r *)' ( 2 ) 

i=l 

where and are the spin components for replicas 
(1) and (2), respectively. It has been proposed that there 
is no SG transition in isotropic Hcisenbcrg SG systems 
but that, instead, the freezing is in the chiral sector—. 
Latest simulations^— suggest that both chiral glass 
(CG) and SG transitions occur at finite temperature, but 
there is no consensus whether the CG critical tempera- 
ture (Tcg) is higher or equal to that of the SG (?sg)- 
Since T-L has isotropic Heisenberg spins, there is no obvi- 
ous reason to exclude the possibility of a CG transition. 
To monitor CG correlations, we consider two chirality 
parameters. The first one is defined along bonds, 

1 N 

9c G1 (k)^— £ £ ^^exp^k-r,), (3) 

% ~ 1 <5=ei,i;«i,2;£i,3 

where K Vi i = S {l ~ s) ■ (Si x S^_^). The e^i, e ij2 and e i>3 
are vectors pointing from site i, to its three nearest neigh- 
bor sites in the same tetrahedron; (i, S), i, (i, —S) are the 
indices for the three sites lying along the direction 5 (see 
Appendix). The second one is defined on the triangular 
faces of individual tetrahedra, 

27V 

<?CG2(k) = — £ k^Ik { ^1 exp(ik ■ r u ), (4) 



where k 2;w = S u . a ■ (S u , b x S WlC ); (w, a), (ui, b), (w, c) 
are the indices for the three sites of the triangular face 
lo and i u = (f U)0 + f w . fa + t UiC )/3. The correspond- 
ing susceptibilities are obtained from the order param- 
eters, XSG (k) = NJ2^ xy J{\ q ^(k)\ 2 )}; XcG1 (k) = 
3N[(\q CG1 (k)\ 2 )]/x?, XcG2 (k) = 2N[(\q CG2 (k)\ 2 )]/^ i , 
where (...) and [...] denote the thermal average and 
disorder average, respectively. As the local chiral- 
ity variables are not fixed to be unity as is the 
case of the spin variable s 19 ! 20 , the CGI susceptibil- 
ity is normalized by xT 4 = (xl 2 ) 2 , where xi 2 = 
W E*LiEj =eM ,i i)2l?i ,3 i( K l-i,s)} ' and > similarly, the CG2 
susceptibility is normalized by xi 4 = (~X2 2 ) 2 , where 
X2 2 = 2lE!!i[( K 2;J] (see Appendix). 

Assuming that the susceptibilities follow an Ornstcin- 
Zernikc fornJd, the correlation lengths £sg and £cgi,cg2 

1 /2 

can be determined via f fl (L) = 2ain( | k | /2) (f^g - l) , 
where g = SG, CGI or CG2 and k = 2-kx/L is one of the 
smallest wave vectors for system size LlZr—. The £ g 's 
divided by L should be scale invariant at their respective 
critical point. The crossing of £, g /L is therefore a sen- 
sitive criterion to test for a glass transition. The corre- 
lation lengths and susceptibilities should finite-size scale 
as £ g {L)/L = X[(T — T g )L 1 / v '>} and Xff (T,L)Z>- 2 = 
Y[(T - Tg)^/""}, respectively. To check that thermody- 
namic equilibrium was reached, we verified that £ g be- 
comes independent of simulation time for the largest sys- 
tem size and lowest temperature considered. 

We first present the SG and CG correlation lengths in 
Fig. [T] It is fairly clear from these results that the SG 
and CG correlation lengths for different system sizes tend 
to cross in a narrow range of temperatures compatible 
with a nonzero critical temperature for both SG and CG 
(T g « 0.016). 

We then employ a scaling scheme which assumes that 
the correlation lengths finite-size scale as the scaling func- 
tion X given above and fit the data in the temperature 
range (0.012 < T < 0.020) by parametrizing X as poly- 
nomials X(z) = J2 m =o.i,...,5 c ™( z ~ z o) m - The merit 

function A, A = Emc data[ X ( z )/fe,i/ L ) - is mini- 
mized numerically to obtain the coefficients c m , Zq, crit- 
ical temperature T g and exponent v g . Figure [5] shows 
£ g (L)/L versus the scaling parameter z = (T — Tg)!, 1 /" 11 
(T g and v g are listed in the caption of Fig. [2]) The 
scaling exponents determined both for the SG and CG 
are far from those of the 3D EA Ising model obtained 
from correlation length scaling (^sg ~ 2.44) 15 i 16 , but 
roughly comparable (within ~ 20%) to those of the 3D 
EA Heisenberg modeU£ ~ 21 i 23 . 

In the above scaling analyses, we assumed that there 
is a common crossing point for all system sizes. Realisti- 
cally, the critical temperature obtained this way for fair 
system sizes (N < 8192) should represent an upper bound 
for the true critical temperature in the thermodynamic 
limit. For example, in the latest simulations of the 3D 
EA Heisenberg model, it was found that scaling correc- 
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FIG. 1: (Color online). Correlation lengths for spin glass 
(top) and chiral glass (middle and bottom). The insets show 
the details close to the crossing points. 



tions are large and that the £ s crossings arc pushed to 
lower, albeit non-zero, temperatures as the system size 
increases^—. We show in Fig. |3] the evolution of the 
correlation lengths crossing temperatures T* for different 
pairs of system sizes (L\, L-i) as a function of the inverse 
of their average size given by 1/ L avc = 2/ (L 1 +L 2 ) . These 
show that < T SG < T CG1 w T CG2 as L avc -> oo. 

We use the same procedure as for the correlation length 
data collapse to fit the susceptibilities into the scaling 



0.6", 

\ 

, 0.5- 

5 0.4- 


^0.3- 
0.2- 



"L=4 
L=8 




-0.02 0.02 0.04 0.06 0.08 0.1 
IN 

(T-T $G )L^ 



o.o L 



0.05 

lAi 

(T-T SG )L SG 



o.i 



0.5- \ 



^ 0.4- 
5,0.3- 

5 

u 02- 
0.1- 
0-_ 

0.5- 
^0.4- 
3" 0.3- 



-L=4 
°L=6 



^ 10 - 

N 





0.01 0.01 0.02 0.03 0.04 0.0! 
(T-T CGl )L 1A ™ 



0.02 0.04 
(T-T CGl )L 1/v ™ 



*L=4 
•L=6 
L=8 



80.2- 
0.1- 



o, 




'0.02 



0.02 0.04 




0.02 0.04 0.06 OX 

(r-r CG2 )z> 



FIG. 2: (Color online). Left column: Correlation lengths 
for spin glass (top) and chiral glass (middle and bottom) 
vs (T — Tg)L 1 ' 1 ' 9 . Critical temperatures and scaling expo- 
nents are T SG = 0.0157, ^sg = 1.037 for SG; T CG i = 0.0153, 
^cgi = 1-585 and T CG2 = 0.0161, u CG 2 = 1.408 for CG, ob- 
tained by fitting the data to the scaling function X. Right 
column: X g(T, L)L V ^ 2 for spin glass (top) and chiral glass 
(middle and bottom) vs (T - T g )L 1/Vg . Critical tempera- 
tures and scaling exponents are Tsg = 0.0126, j^sg = 1.024, 
?7sg = -0.292 for SG; T CG i = 0.0134, v C Gi = 1.439, r, CG i = 
0.621 and T CG 2 = 0.0139, i/ C G2 = 1.155, 7? C G2 = 0.698 for 
CG, obtained by fitting the data to the scaling function Y. 
The numbers listed are the actual values used in these scaling 
plots. Corrections to scaling are expected to be sizable, there- 
fore no error bar is provided, while the exponents lie within 
the range estimated for the 3D EA Heisenberg model. 



function, Y, in order to determine T 9 , v g and r\ g . Figure 
[5] shows Xg(T, L)L Vg ~ 2 versus the scaling parameter z = 
(T - T g )L 1 l v o. (The T g , v g , and r/g are listed in the 
caption of Fig. [2]) The rj g and v g values are again fairly 
comparable with those obtained in the latest studies of 
the 3D EA Heisenberg modet^ -21 ' 23 , providing evidence 
for a common SG universality class for the BGHM model 
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FIG. 3: (Color online). Left: Correlation length crossing tem- 
peratures for SG, CGI and CG2 for different pairs of system 
sizes (Li, L2). The vertical axes are the crossing temperatures 
and the horizontal axis is the inverse of the average system 
sizes given by 1/L ave = 2/(Li + L2). The lines are guides for 
linear extrapolations to the limit L avc — ► oo, that is assum- 
ing Tg(Li, L2) — T g oc L^ve with = 1. Right: The average 
tetrahedra moment, <£>, vs temperature, T. 



and the 3D EA Heisenberg model. 

A recent study aimed at describing the SG in H as- 
sumes that the power-law correlation of the CSL^ is 
maintained despite the random disorder and thus the 
spins can be thought of as interacting via an effective 
projected interaction matrix of a long range "dipolar" 
for m 11 ' 12 as a consequence of the zero net magnetic mo- 
ment ($t=0) condition on each tetrahedron. To inves- 
tigate this description, we calculate the average tetra- 
' 12 1 1 Ei=i 4 s (t,o I ' 



hedra moment, $ 



where the 



outer sum is over all tetrahedra, the inner sum is over the 
four spins in each tetrahedron, and -/Vt etr ahedron = N/2 is 
the total number of tetrahedra. We show $ as a func- 
tion of temperature in Fig. [3] First, we find that it 
changes very little with system size, as it is not a criti- 
cal quantity. Second, it decreases with decreasing tem- 
perature. Most importantly, near the crossing tempera- 
tures (T g w 0.016), $ is finite and of the order of W/J 
which implies the existence of "defect" tetrahedra with 
$t 7^ and, consequently, the destruction of infinite- 
range power-law correlations^. This is further supported 
by the reasonably good data collapse for the £ g corre- 
lation lengths extracted from an Ornstein-Zernike form, 
which would likely not be correct if there were a CSL 
phase with extended power-law correlations intervening 
between the paramagnetic phase and the SG phase. 

For Y2M02O7, the experimentally determined Curie- 
Weiss and SG temperatures are O^,' « —200 K 2 ^ and 
TgQ P ' « 22.5 K^, respectively. The nearest neighbor co- 
ordination number on the pyrochlore lattice is z = 6 and 
the spin S of magnetic Mo 4+ is S = 1. Therefore, with 
J ~ J^ S 2 and = J™ ic zS(S + l)/3 « 200K, we 

get J mic - ps 50 K and the ratio T|^/(J mic -S 2 ) « 0.45. 
This is much higher than that obtained for the BGHM 
model for which we found above Tsg/Jo ~ 0.01 — 0.02 
for W/Jq = 0.1. This suggests that the glass transition 



of Y2M02O7 is not due to weak random disorder as in H, 
but rather to very strong effective disorder. One plausible 
scenario is that perturbations beyond nearest-neighbor 
Heisenberg exchange Jo disrupt the perfect degeneracy of 
the CSL phase and induce short range AFM order above 
T 9 , as observed in a neutron scattering study 3 -. If the 
growth of AFM order is forestalled due to some form of 
random disorder— ~—, the SG behavior of Y2M02O7 should 
likely be described in terms of a "cluster-glass" modcli. 

In conclusion, our MC simulations of the BGHM model 
of Eq. (1) provide compelling evidence for a thermody- 
namic SG phase induced by weak random disorder in a 
classical spin liquid of a highly frustrated system. From 
our work, it appears very likely that the SG behavior 
in Y2M02O7 is not due to weak and dense random dis- 
order but rather via an effective strong disorder whose 
microscopic origin requires further investigation. 

This work was funded by the NSERC of Canada, the 
Canada Research Chair Program (M.G., Tier 1) and 
SHARCNET. We thank H. Kawamura for encouraging 
us in studying this problem and for his useful comments 
on our manuscript, and J. Chalkcr, P. Holdsworth, P. 
McClarty and P. Stasiak for useful discussions. 



Appendix A: Spin and Chiral Glass Susceptibilities 

This appendix discusses the details of the definitions 
of the spin glass (SG) and two different chiral glass (CG) 
overlap parameters - CGI defined along the bonds and 
CG2 defined on the triangular faces of the tetrahedra that 
from a pyrochlore lattice. In addition, a table for critical 
exponents obtained from previous Monte Carlo studies 
of the three-dimensional (3D) Edwards- Anderson (EA) 
Heisenberg SG model is provided for comparison with the 
critical exponents obtained for the BGHM model studied 
in this paper. 

The spin glass (SG) overlap is defined as the overlap 
between two thermal replicas with the same realization 
of random couplings {Jfj S }, 



N 



^G(k)^VS^exp(zk. 



■ =1 



(Al) 



where and sf*^ are the spin components of the two 
replicas. For Ising spins, this is the usual parameter used 
to monitor the spin freezing. The situation is, however, 
more complicated for Heisenberg spins^. 

It has been proposed that there is no SG transition in 
isotropic Heisenberg SG systems but, instead, that the 
freezing is in the chiral sector—. There have been many 
investigations on the spin-chirality coupling/decoupling 
scenario for the 3D EA Heisenberg model ^ 22 ' 24 " — , but 
there is so far no consensus ^ 2 1 ' 30 . While recent studies 
tend to agree that both the critical temperature for SG 
(Tsg) and chiral glass (Xcg) are non-zero 1 ^—, whether 
T CG = T so ^- or T CG > T SG ^. is still under active 
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debate. Since the model we study possesses spins with 
isotropic Hcisenberg exchange, there is no a priori reason 
to rule out the possibility of a chiral glass (CG) transi- 
tion. 

To monitor the CG we probe two different chirality 
variables. The first one is defined along the bonds. This 
is a generalization of the definition of chirality variables 
for the 3D EA Hcisenberg model on a simple cubic lattice. 
Similarly to the simple cubic lattice, three are three axes 
passing through each site. As three spins are needed 
to define the chirality, the natural choice is to pick a 
spin and its two nearest neighbors spins along one of the 
axes to define a chirality variables. As there are three 
axes passing thorough each site, therefore there are three 
chirality variables for each site. We denote this definition 
of chirality overlap as CGI. 

1 N 

where K 1;i j = ■ (S l x S^_^). The e^i, e ij2 and 

£i.3 arc vectors pointing from site i, to its three near- 
est neighbor sites in the same tetrahedron, see Fig. |U 
(i, 6), i, (i, — o) are the indices for the three sites lying 
along S. The normalization factor 1/37V is introduced to 
account for the 3N chirality variables for this definition 
because of the three chirality variables at each site. We 
note that a slightly different definition, which treats the 
CG overlap as a three components object has been used 
in a study of 3D EA Heisenberg model^i. 




FIG. 4: (Color online). The figure illustrates the vectors for 
defining the chirality at one of the sites. The blue color arrows 
show the direction of the vectors ej,i, e;,2 and ei,3 for the site 
i (in the middle) shared by two tetrahedra. 

The second chirality overlap parameter is defined on 
the triangular faces of the tetrahedra. The pyrochlore lat- 
tice is composed of corner-sharing tetrahedra, and each 
tetrahedron has 4 triangular faces. Therefore, another 
natural choice is to define the chirality variables on each 
face of the tetrahedra. In a lattice with N sites, there are 
N/2 tetrahedra, and each tetrahedron has 4 faces. There- 
fore, there are in total 2N chirality variables within this 
definition. We denote this definition of chirality as CG2, 



with: 

<?CG2(k) = — £ exp(ik ■ i u ), (A3) 

where k 2;w = S U) a ■ (S u ,b x Sk,c); (w, a), (w, 6), (w, c) are 
the indices for the three sites of the triangular face ui 
and f w = (fa;, a + ?a;, 6 + r WjC )/3. The normalization fac- 
tor I /27V is introduced to account for the 2N chirality 
variables for this definition since there are 2N triangular 
faces. 

The corresponding susceptibilities are obtained from 
the above order parameters, 

XsG (k) = N Kl9sG(k)| 2 )]; (A4) 

fj,,v—x,y,z 

XcGi(k) = 3iV[(| gcG1 (k)| 2 )]/xr 4 , (A5) 
XcG2 (k) = 2/v[(| (7 cG2(k)| 2 )]/x2- 4 , (A6) 

where (...) and [...] denote the thermal average and dis- 
order average respectively. As the local chirality vari- 
ables are not fixed to be unity in contrast of the spin 
variable s 19 ' 20 ' 35 , the CGI susceptibility is normalized by 
XT 4 = (Xl 2 ) 2 > where 

1 N 

xr 2 -^E_ E [<«?„,«>]. ( A7 ) 

and the CG2 susceptibility is normalized by xi 4 = 
(%2 2 ) 2 , where 

2N 

^^EtHJ]- (A8) 

ui=X 

The exponents obtained for the BGHM model are 
fairly comparable with that of the 3D EA Heisenberg 
spin glass model 3 ^. For comparison of our results with 
the 3D EA Heisenberg model, we compile a selection of 
critical temperatures and critical exponents in Table ITT] 
For the 3D EA Ising model, see Table I in Ref. 
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reference 


randomness type 




"SG 


VSG 




"CG 


7?CG 


Kawamura^ 


Gaussian 


NA 


NA 


NA 


0.157 ±0.01 


NA 


NA 


(1998) 
















Hukushima and Kawamura^ 


Gaussian 


NA 


NA 


NA 


0.160 ±0.005 


1.2 


0.8 


(2000) 
















Endoh, et aZ.— 


Bimodal 


0.19 ±0.02 


NA 


NA 


NA 


NA 


NA 


(2001) 
















tT"t i 'ju 

Matsubara, et at.— 


Bimodal 


0.18 


NA 


NA 


NA 


NA 


NA 


(2001) 
















Nakamura and Endoh^ 


Bimodal 


u -^ 1 -0.03 


1.1 ±0.2 


0.27* 


u -^ z -0.04 


NA 


NA 


(2002) 












Lee and Youngi^ 


Gaussian 


0.16 ±0.02 


1.1 ±0.2 


NA 


Tcg = Tsg 


1.3 ±0.3 


NA 


(2003) 
















Nakamura, et a/.— 


Bimodal 


0.20 ±0.02 


0.8 ±0.2 


-0.375* 


NA 


NA 


NA 


(2003) 
















Hukushima and Kawamura^ 


Bimodal 


NA 


NA 


NA 


0.19 ±0.01 


1.3 ±0.2 


0.8 ±0.2 


(2005) 
















Viet and Kawamura^ 


Gaussian 


n ioc+U.UUB 
U.±ZO_q 012 


NA 


< -0.30 


0.143 ± 0.003 


1.4 ±0.2 


0.6 ±0.2 


(2009) 














Fernandez, et ali^ 


Gaussian 


120 +u ' uiu 
u.±zu_ 004 


1.49 ±0.13* 


-0.19 ±0.02* 


Tcg = Tsg 


1.30 ±0.08* 


0.56 ±0.04* 


(2009) 














This work (£,/L scaling) 


pyrochlore 
±0.1Jo uniform 


0.0157 


1.037 


NA 


0.0153(CG1) 
0.0161(CG2) 


1.585(CG1) 
1.408(CG2) 


NA 


This work (x scaling) 


pyrochlore 
±0.1Jo uniform 


0.0126 


1.024 


-0.292 


0.0134(CG1) 
0.0139(CG2) 


1.439(CG1) 
1.155(CG2) 


0.621(CG1) 
0.698(CG2) 



TABLE II: Selection of critical temperatures and exponents of the three-dimensional Edwards- Anderson Heisenberg model for 
the simple cubic lattice. For the Edwards- Anderson Ising model see Table I in Ref. [l5|]. NA is a shorthand for not available. 
The last two rows are the estimates from this work. The variance of the random coupling distribution is 1 for all the studies 
on the three-dimensional Edwards- Anderson Heisenberg model listed in this table, and the variance of the model we study is 
1/3 x 10~ 2 . (The variance is defined as f(x — x) 2 P(x)dx, where P(x) is the distribution function and x — f xP(x)dx.) The 
entries with are not quoted in the original paper but are estimated via the relation 7 = (2 — r\)v. The entries with "*" are 
from the quotient metho d 41 ' 42 for scaling analysis of L — 24 and L = 48 systems on a simple cubic lattice. 
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